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ABSTRACT 



^> , We present moderate- resolution spectroscopy of globular clusters around the Sa/Sb spiral 

■ galaxy M81 (NGC 3031). Sixteen candidate clusters were observed with the Low Resolution 

5^ , Imaging Spectrograph on the Keck I telescope. All are confirmed as bona fide globular clusters, 

although one of the clusters appears to have been undergoing a transient event during our ob- 
servations. In general, the M81 globular cluster system is found to be very similar to the Milky 
Way and M31 systems, both chemically and kinematically. A kinematic analysis of the velocities 
of 44 M81 globular clusters, (the 16 presented here and 28 from previous work) strongly suggests 
that the red, metal-rich clusters are rotating in the same sense as the gas in the disk of M81. The 
blue, metal-poor clusters have halo-like kinematics, showing no evidence for rotation. The kine- 
matics of clusters whose projected galactocentric radii lie between 4 and 8 kpc suggest that they 
are rotating much more than those with projected galactocentric radii outside these bounds. We 
suggest that these rotating, intermediate-distance clusters are analogous to the kinematic sub- 
population in the metal-rich, disk globular clusters observed in the Milky Way and we present 
evidence for the existence of a similar sub-population in the metal- rich clusters of M31. With 
one exception, all of the M81 clusters in our sample have ages that are consistent with Milky 
Way and M31 globular clusters. One cluster may be as young as a few Gyrs. The correlations 
between absorption- line indices established for Milky Way and M31 globular clusters also hold 
in the M81 cluster system, at least at the upper end of the metallicity distribution (which our 
sample probes). On the whole, the mean metallicity of the M81 globular cluster system is similar 
to the metallicity of the Milky Way and M31 globular cluster systems. To within a factor of two, 
the projected mass of M81 is similar to the masses of the Milky Way and M31. Its mass profile 
indicates the presence of a dark matter halo. 



Subject headings: galaxies: spiral, galaxies: individual (M81,NGC3031,M31), globular clusters: general 



1. Introduction 

Globular clusters arc ideal probes of the dynam- 
ical properties and chemical histories of their host 
galaxies. They are bright enough to be observed 
at large distances, old enough to have been "wit- 
nesses" to the processes by which galaxies form 
and evolve, and are found in all types of galax- 
ies. For all but the nearest galaxies, globular clus- 
ters are the sole providers of information about 
the stellar population in galaxy halos. 

Chemically and kinematically distinct sub- 
systems have been established in the globular clus- 
ter systems of both the Milky Way and our nearest 
spiral neighbor, M31. The presence of such sub- 
systems places significant constraints on spiral 
galaxy formation theories and may also bear on 
understanding the formation and evolution of el- 
liptical galaxies. It has been suggested (Schweizer 
1986) that some elliptical galaxies may form via 
the merging of spirals. In a gaseous merger, the re- 
sulting globular cluster system would be expected 
to contain any globular clusters originally asso- 
ciated with the merger participants and perhaps 
new clusters formed during the merger (Ashman 
& Zepf 1992). While study of globular cluster sys- 
tems around ellipticals is important for validating 
or refuting the merger hypothesis, interpreting 
such studies requires understanding the natural 
variations in the observed properties of globular 
cluster systems around spirals. Key properties 
such as age and metallicity are most reliably con- 
strained through spectroscopy, while kinematic 
and individual element abundance information are 
obtained exclusively so. The advent of 10-meter 
class telescopes and multi-object spectrographs 
has made progress on this front feasible. 

The Milky Way globular cluster system is natu- 
rally the benchmark by which the systems of other 
spirals are measured. One of its most striking fea- 
tures is its bimodal metallicity distribution, with 
peaks at [Fe/H] = -1.6 and -0.5 (Zinn 1985; Ar- 
mandroff & Zinn 1988; Armandroff 1989). Fur- 
thermore, these metallicity sub-populations are 
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spatially and kinematically distinct (Zinn 1985; 
Hesser, Shawl & Meyer 1986; Armandroff & Zinn 
1988; Armandroff 1989). The metal-rich clusters 
comprise a centrally-concentrated, flattened, ro- 
tating system with a low velocity dispersion, re- 
ferred to by Zinn (1985) as the disk population. 
Armandroff (1989) associated these disk clusters 
with the Milky Way's thick-disk component, based 
on similar spatial, kinematic and chemical prop- 
erties of the two systems. More recently, it has 
been reported that the properties of the metal-rich 
clusters lying within 3kpc of the Galactic center 
have more in common with the bulge than with the 
thick disk (Minniti 1995; Cote 1999). Burkert & 
Smith (1997) further explored sub-structures in 
the metal-rich globular cluster population, identi- 
fying a group belonging to the galactic bar and one 
comprising a ring of rapidly-rotating clusters with 
a low velocity dispersion and galactocentric radii 
between 4 and 6 kpc (these are the thick-disk clus- 
ters of Armandroff (1989)). To first order, Zinn 
(1985) found the properties of the Milky Way's 
metal-poor globular clusters to be well matched 
with those of the halo field stars: they are spheri- 
cally distributed, rotate slowly and have a large 
line-of-sight velocity dispersion. However, Zinn 
(1993) reported substructure in the Milky Way's 
halo clusters as well. 

The M31 globular cluster system was the first 
system around another spiral to be studied in de- 
tail. A number of imaging programs resulted in 
several catalogues and refined lists of globular clus- 
ter candidates (Crampton et al. 1985; Battistini, 
et al. 1980, 1987; Battistini et al. 1993; Reed, Har- 
ris & Harris 1994). Insight into the chemical and 
kinematic nature of M31's globular cluster sys- 
tem resulted from the optical spectroscopy stud- 
ies of Huchra, Stauffer & Van Speybroeck (1982), 
Kent, Huchra & Stauffer (1989), Huchra, Brodie 
&Kcnt (1991) and Brodie & Huchra (1991). All 
came to the same conclusion: the globular cluster 
system of M31 is, for the most part, very similar 
to the Milky Way system. Huchra, Brodie & Kent 
(1991) derived metallicities for 149 clusters and 
found them to be consistent with the metallicity 
distribution of the Milky Way globular cluster sys- 
tem, but they were unable to distinguish any obvi- 
ous bimodality. They noted that large uncertain- 
ties in their individual metallicity measurements 
could be obscuring its presence. Ashman & Bird 
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(1993) subsequently performed a statistical analy- 
sis on the Huchra, Brodie & Kent (1991) data and 
found quantitative evidence for a bimodal metal- 
licity distribution at the 98.4% confidence level. 
New optical and infrared color distributions by 
Barmby et al. (2000) definitively support the pres- 
ence of bimodality. 

Based on velocities of 60 clusters, Huchra, 
Stauffer & Van Speybroeck (1982) detected a 
modest, but significant rotation for the entire sys- 
tem. Subsequent studies have yielded evidence for 
kinematic sub-systems as well. Huchra, Brodie & 
Kent (1991) reported that the metal-rich clusters 
in the inner few kpc of M31 are rapidly rotating, 
while at larger radii they observed little rotation 
in the M31 globular clusters, regardless of cluster 
metallicity. Barmby et al. (2000) recently re- 
peated the kinematic analysis of Huchra, Brodie 
& Kent (1991) with an additional 74 clusters and 
derived the same result. 

As similar to the Milky Way system as M31's 
globular cluster system appears, there do seem to 
be some genuine differences. The M31 clusters 
contain, on average, a higher abundance of nitro- 
gen than Milky Way globular clusters. This result 
was first reported as a CN enhancement in the 
spectroscopic studies of Burstein et al. (1984) and 
Brodie & Huchra (1991) and more recently in the 
ultraviolet spectroscopy of Ponder et al. (1998), 
who established it as a nitrogen overabundance. 
Brodie & Huchra (1991) also reported evidence 
of calcium enhancement in the M31 cluster spec- 
tra. In addition, the M31 system contains a hand- 
ful globular clusters which appear to be "young" , 
based on the strength of their Balmer-line ab- 
sorption (Sargent et al. 1977; Cowley & Burstein 
1988; Elson & Walterbos 1988; Brodie & Huchra 
1991; Barmby et al. 2000). No such clusters have 
been observed in the Milky Way. 

Globular cluster systems around a few other 
spirals have been examined to a more limited ex- 
tent. Photometric (Christian & Schommer 1982; 
Cohen, Persson & Searle 1984; Sharov & Lyu- 
tyj 1984; Christian & Schommer 1988; Chandar 
et al. 1999a,b,c) and spectroscopic (Schommer et 
al. 1991; Brodie & Huchra 1991) studies of M33, 
the late-type, low-luminosity, Local-Group spiral, 
have revealed that it contains a number of young 
and intermediate-age compact star clusters. In ad- 
dition, it contains about 25 typical old globular 



clusters in a halo distribution with little net rota- 
tion. 

Beyond the Local Group, M104 and M81 are 
the only spiral galaxies whose globular cluster sys- 
tems have been studied spectroscopically. Veloci- 
ties for 34 Ml 04 globular clusters were obtained 
by Bridges et al. (1997), who found marginal 
evidence for rotation. The quality of their data 
was too poor to measure metallicities for individ- 
ual clusters. They produced a rough estimate of 
the mean metallicity of the system from a compos- 
ite spectrum of all 34 objects, finding ([Fe/H]) = 
— 0.7±0.3. This value is substantially higher than 
the metallicity of the Milky Way and M31 globu- 
lar cluster systems, but Ml 04 is extremely bulge- 
dominated. It is also significantly more luminous 
than the Milky Way or M31, and its mean globu- 
lar cluster metallicity is similar to the values found 
for globular clusters around elliptical and giant el- 
liptical galaxies with a luminosity similar to that 
of M104 (see Table 4 in Bridges et al. 1994.) 

At a distance of 3.7 Mpc (Preedman et al. 1994) 
M81 is at roughly half the distance of M104. It is 
an Sa/Sb-type spiral, very similar to M31, and 
roughly as massive as the Milky Way. Its proxim- 
ity and morphology make it a good candidate for 
reaching beyond Local Group toward a detailed, 
high signal-to-noise spectroscopic study of a spi- 
ral galaxy globular cluster system for comparison 
to the Milky Way and M31 systems. In their 
study of extragalactic globular clusters, Brodie & 
Huchra (1991) derived spectroscopic metallicities 
for eight clusters in M81 and found a sample mean 
of [Fe/H] = -1.46 ± 0.31. Perelmuter & Racine 
(1995) conducted a photometric study of the M81 
globular cluster system and compiled a catalogue 
of globular cluster candidates, using proper mo- 
tion to minimize contamination from foreground 
stars. In a spectroscopic follow-up, Perelmuter, 
Brodie & Huchra (1995) obtained low signal-to- 
noise spectra of 82 candidates, 25 of which were 
established as bona fide members of the M81 glob- 
ular cluster system. The generally poor signal-to- 
noise of their data resulted in rather large uncer- 
tainties in the metallicities of individual clusters. 
They derived a more precise mean metallicity of 
their sample from a composite spectrum, finding 
([Fe/H]) = -1.48±0.19. Their analysis of the clus- 
ter velocities showed a hint of rotation in the M81 
globular cluster system. 
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In this work wc extend the spectroscopic stud- 
ies of Brodie & Huchra (1991) and Perelmuter, 
Brodie & Huchra (1995) with significantly higher 
signal-to-noise spectra of an additional 16 mem- 
bers of the M81 globular cluster system. Section 
2 describes our observations and data reduction 
procedure. A discussion of the abimdance ratios, 
mean metallicities and ages of the M81 clusters 
and a comparison of these properties to M31 and 
the Milky Way globular clusters can be found in 
Section 3. In Section 4, we combine our veloc- 
ity data with the velocities of additional clusters 
from the work of Perelmuter, Brodie & Huchra 
(1995) and Brodie & Huchra (1991), and discuss 
the kinematics of the resulting sample. Section 
5 contains a summary of our results and a brief 
discussion of their implications for using globular 
cluster systems as probes for studying galaxy for- 
mation and evolution. 

2. Observations and Data Reduction 
2.1. Observations 

The targets were selected from the candidate 
list of Perelmuter & Racine (1995). The specific 
candidates observed were selected on the basis of 
brightness and the desire to simultaneously ob- 
serve the maximum possible number of candidates 
with the spectrograph in multi-object mode. Spec- 
tra of 16 globular cluster candidates were obtained 
on 1994 November 30 using the Low Resolution 
Imaging Spectrograph (LRIS) (Oke et al. 1995) 
on the Keck I 10-m telescope in ~ 0."8 seeing. 
The candidates were obscirvcid through a single slit 
mask. Weather and technical problems prevented 
observations through two additional masks, so our 
sample is spatially biased toward one side of the 
galaxy. The 600 1/mm grating, blazed at 7500A, 
and slitlets of 1" width resulted in a dispersion 
of 1.25 A/pixel and a spectral resolution of 5-6A. 
Because the horizontal positions of the slitlets var- 
ied across the slit mask, the spectral range falling 
on the detector was different for each slitlet. The 
shortest wavelength observed was 3700A and the 
longest was 7420A. Each candidate spectrum cov- 
ers 2560 A between those two extremes. The re- 
gion of common overlap for all candidates is 4680- 
6260A. We obtained 5 x 900 s integrations for a 
total integration of 4500 s. A complete list of 
the observed cluster candidates' positions, magni- 



tudes and colors is provided in Table 1. Figure 1 
shows the positions of the candidates with respect 
to the galaxy center, along with the positions of 
established chistcrs from the work of Perelmuter, 
Brodie & Huchra (1995) and Brodie & Huchra 
(1991). 

2.2. Data Reduction 

These data were reduced using the following 
procedure. A mean bias, assumed to be the mean 
level of the CCD overscan region, was subtracted 
from both science and calibration images. In addi- 
tion, the small-scale bias structure was determined 
by combining multiple one-second dark exposures. 
This bias structure was subtracted from all other 
science and calibration images. After the images 
were cleaned of cosmic rays, a map of the small 
scale variations in detector sensitivity and grat- 
ing efficiency was removed by fitting a fifth- or 
sixth-order polynomial to a halogen lamp spec- 
trum taken through each slitlet in the mask and 
dividing the lamp spectra by the polynomial. The 
resulting normalized images were divided into the 
science images to remove the small scale varia- 
tions. 

Each slitlet was subsequently processed using 
the IRAF^ routines of the noao .twodspec pack- 
age. Ident, reident and f itcoords were run on 
spectra of Hg, Ne, Ar and Kr lamps to determine 
the wavelength solution and then transform was 
used to apply the appropriate solution to the two- 
dimensional spectra of the candidate objects. The 
object spectra were extracted using the IRAF rou- 
tines of the noao . apextract package, employing 
the optimal spectral extraction scheme detailed by 
Horne (1986). The polynomial used to approxi- 
mate the background was typically a linear fit. For 
a few candidates, however, proximity to end of the 
slitlet required approximating the background as 
a constant. 

The extracted candidate spectra were flux- 
calibrated using a long slit spectrum of the stan- 
dard star EG145 taken on the same night as the 
candidate observations. Since its location in the 
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dispersion direction on the mask determines what 
wavelength region falls onto the detector for an in- 
dividual slitlet, the wavelength coverage for some 
candidate cluster spectra is not coincident with 
the coverage for the flux standard. The result is 
that the continuum shape is only accurately de- 
termined in the spectral region of 3900A-6500A. 
None of the absorption features of interest lie be- 
yond 6500A, while two spectral features of inter- 
est (sec Section 3) require measurements below 
3900A. The final spectra of all 16 candidates are 
shown in Figure 2. They have been corrected 
for foreground reddening in the direction of M81 
using the IRAF task noao . onedspec . deredden 
with E(B-v) =0.1 and Ry = 3.1 (Burstein & Heiles 
1984), then shifted to the rest frame. 

M81's low systemic heliocentric velocity (—34 km s 
precludes the establishment of globular clus- 
ter candidates as bona fide M81 globular clus- 
ters on the basis of radial velocity alone. Perel- 
muter, Brodie & Huchra (1995) circumvented this 
problem by comparing each candidate's color to 
the relative strengths of the calcium K and B.6 ab- 
sorption lines in its spectrum. Unfortunately, the 
low throughput of LRIS below 4000A and the non- 
coverage of that part of the spectrum for many of 
our candidates prevented us from doing likewise. 
However, since our field is within 5' of the center 
of M81 where the cluster surface density is high, it 
is very likely that candidates whose spectra con- 
tain typical globular cluster features are members 
of the M81 globular cluster system. This asser- 
tion is supported by the fact that the Bahcall & 
Soneira (1980) models predict fewer than 2 galac- 
tic stars in our 48 arcmin^ field (the LRIS field is 
6'x 8') in the direction of M81 and in the color 
and magnitude range of our candidates. 

Fifteen of the 16 candidates observed have spec- 
tra that are typical of globular clusters and we 
henceforth treat them as such. Only one, object 
12, appears at first to be ruled out as a globular 
cluster on the basis of its spectrum (see Figure 2). 
The sharp upturn at the blue end of the spectrum 
and what appears to be broad-line Balmer emis- 
sion are decidedly "un-globular" . Still, there are 
genuine globular cluster features in the spectrum 
of this object, occurring at wavelengths that are 
consistent with it belonging to M81. Curiously, 
the photometric study of Pcrelmuter & Racine 
(1995) cataloged a B — V value of 1.0 for this ob- 



ject, which is well within the normal color range 
for a globular cluster. The color of an object with 
a spectrum such as object 12 would certainly be 
much bluer, and likely would not have satisfied the 
Perelmuter & Racine (1995) color criteria for in- 
clusion in their catalogue of candidates. Thus we 
suspect that a transient phenomenon was occur- 
ring in the cluster at the time of our observation. 
Preliminary reductions of a spectrum of object 
12 recently obtained with the Echellettc Spectro- 
graph and Imager (ESI) on the Keck II telescope 
has confirmed our suspicion. Object 12 no longer 
shows an upturn at the blue end of the spectrum 
and no broad-line emission features are present. 
Thus we consider object 12 to be a genuine M81 
globular cluster and include it in our kinematic 
1) analysis of the M81 globular cluster system (Sec- 
tion 4). However, because the continuum shape 
is pivotal for accurate measurements of absorp- 
tion line strengths and hence metallicity, we ex- 
clude object 12 from our chemical abundance and 
metallicity analysis in Section 3. The exact nature 
of the phenomenon occurring during the original 
observation of object 12 will be the subject of a 
forthcoming paper. 

2.3. Measurement of Radial Velocities 

A radial velocity was determined for each clus- 
ter by cross-correlating its non-flux-calibrated 
spectrum with three high signal-to-noise tem- 
plates of M31 globular clusters. The actual cross- 
correlation was accomplished using the IRAF rou- 
tine f xcor in the rv package. One of the templates 
(225-280) was obtained on the same run as these 
data and the other two (225-280 and 158-213) were 
obtained on a subsequent run, also with LRIS and 
at the same resolution. The ID numbers and helio- 
centric velocities of the template objects are from 
Huchra, Brodie & Kent (1991). The template 
spectra have the same resolution as the cluster 
spectra, 5-6A, corresponding to a velocity resolu- 
tion of roughly 300 km s^^. The results from the 
three templates agree quite well, the differences 
being a small (^5-10%) fraction of the velocity 
uncertainties returned by f xcor. The heliocentric 
line-of-sight velocities presented in Table 2 are the 
average of the results from the three templates. 
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3. Abundances and Metallicities 
3.1. Metallicities 

We derived each cluster's mean metallicity by 
taking the error-weighted average of nictalhcities 
predicted by the strength of several absorption fea- 
tures. 

The absorption line indices were measured as 
described in Brodie & Huchra (1990), but our 
method for estimating the uncertainties was han- 
dled slightly differently because our background 
subtraction technique was different. 

The extraction routines we used (see section 
2.2) allow for the output of a "variance" spectrum, 
which is actually the standard deviation for each 
pixel in the extracted, background-subtracted 
spectrum, based on the quality of the polyno- 
mial fit to the background and the number of 
rows over which the spectrum was summed before 
extraction. We used this "variance" spectrum to 
estimate the photon uncertainty on the average 
flux per angstrom in a given bandpass by sum- 
ming the "variance" pixels in quadrature over en- 
tire bandpass and then dividing by the bandpass 
width: 



(O-BP) = 




where {(Tbp) is the photon imccrtainty on the 
average flux per angstrom in a generic bandpass, 
apix is the standard deviation in a given pixel (as 
given by the "variance" spectrum), n is the num- 
ber of pixels in the bandpass and A2 and Ai are 
the lower and upper wavelength limits of the band- 
pass, respectively. 

As in Brodie & Huchra (1990), we computed 
the uncertainty in the measured value of an in- 
dex in a particular spectrum, ai, by summing in 
quadrature the photon uncertainties on the flux of 
each bandpass used to compute the index, includ- 
ing the absorption feature bandpass itself and the 
two continuum bandpasses: 

af = +(tI+ (7^2 > 

CI and C2 stand for the two continuum band- 
passes and L stands for the absorption line band- 
pass. 



For each cluster, both the indices themselves 
and the photon errors on the indices were com- 
puted separately for each of our flve integrations. 
The final index values are a weighted average of 
the indices derived from each integration using ^ 
as the weight. The final index photon errors, aj, 
are given by the reduced photon error: 



Table 3 gives the metallicity predicted by each 
spectral feature and the error-weighted mean 
metallicity of each cluster. The uncertainty in the 
metallicities predicted by a given individual index 
were derived by simply propagating the error in 
that index, aj , through the linear relationships be- 
tween the index and metallicity derived in Brodie 
& Huchra (1990). The uncertainty in the error- 
weighted mean metallicity for each cluster is, in 
a sense, an internal error since it is the standard 
of the mean. The individual Brodie & Huchra 
index values may be found in Table 4. Two of 
the features used by Brodie & Huchra (1990) 
to predict metallicity (CNB and Delta) require 
measurements below 3900A. Since, as mentioned 
in Section 2, the continuum shape of the spectra 
is uncertain for regions below 3900A and above 
6500A we measured the index value and predicted 
metallicities for these two features (when the spec- 
tral coverage was siifficicnt) but excluded them 
from our derivation of the clusters' mean metal- 
licities. 

The metallicities of our sample clusters span a 
range of —1.49 < [Fe/H] < —0.21 with an error- 
weighted mean of ([Fe/H]) = -0.91±0.09 and a 
dispersion of 0.29±0.07 dcx. This range is nar- 
rower than that spanned by either the Milky 
Way or M31 globular cluster systems. This may 
be due to a bias in our observed sample. As shown 
in Figure 3 the colors of our sample clusters are 
biased toward the red end of the color distribu- 
tion. The mean projected galactocentric radius of 
our sample is 2.6 kpc. If a cluster population asso- 
ciated with a metal-rich bulge or thick disk exists 
in M81, then our luminosity-selected sample taken 
from a 6' X 8' region centered at 2.6 kpc is likely 
dominated by such clusters, and the mean metal- 
licity of our observed sample is not likely to accu- 
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ratcly represent the mean metallicity of the entire 
M81 cluster system. Indeed, the mean projected 
galactocentric radius of the Perelmuter, Brodie & 
Huchra (1995) sample is significantly larger (9.6 
kpc) and yields a rather lower mean metallicity 
value of [Fe/H] = — 1.48 ± 0.19 as measured from 
a composite spectrum. Brodie & Huchra (1991) 
derived nearly the same result. To within the er- 
rors, this value is consistent with the peak of the 
metal-poor population in the metallicity distribu- 
tions of the Milky Way and M31 cluster systems 
(—1.5 and —1.6, respectively). The mean metal- 
licity of the M81 Keck sample is somewhat lower 
than the mean of the metal-rich peak of M31 and 
Milky Way cluster metallicity distributions. With 
such a small sample size we cannot say whether 
the metal-rich M81 clusters, as a group, have a 
lower metallicity than their counterparts in the 
other systems or if we are simply seeing the inclu- 
sion of a few clusters which belong to a metal-poor 
halo but have small projected galactocentric radii. 

The absence of clusters in our sample with 
[Fc/H] greater than —0.2 is not really surpris- 
ing, given the rarity of such clusters in the Milky 
Way and M31. Only 7 of the 150 M31 clusters 
studied by Huchra, Brodie & Kent (1991) have 
metallicities higher than —0.2, and just three of 
the 134 Milky Way clusters listed in the McMas- 
ter's catalog (Harris 1996) fit this criteria. 

In an attempt to estimate the mean metallicity 
of the entire M81 globular cluster system we de- 
rived the arithmetic mean of all the M81 clusters 
for which spectroscopic metallicities have been es- 
tablished. An error-weighted mean would not be 
useful here because the metallicity errors for the 
Keck sample are miich smaller than the errors in 
the other two studies and the clusters from the 
Keck sample woiild entirely dominate the result. 
In Table 3 of their paper, Perelmuter, Brodie & 
Huchra (1995) indicate that three of the clusters 
in their study overlap with the Brodie & Huchra 
(1991) sample. After comparing the coordinates 
and colors of objects in the two studies we have 
identified all of the objects from Brodie & Huchra 
(1991) as candidates in the Perelmuter & Racine 
(1995) catalogue (see Table 5). There are indeed 
three clusters which are present in both studies 
but one of them was misidentified by Perelmuter, 
Brodie & Huchra (1995). Using the nomencla- 
ture of their Table 4, it is cluster Id50415 that 



should be matched with G19, not Is50401. At any 
rate, we adopt the error-weighted average of the 
metallicities derived in the two studies for these 
three clusters. The arithmetic mean metallicity of 
the resulting sample of 44 clusters is ([Fe/H]) = 
— 1.25 ± 0.13 where the error is the standard error 
of the mean. The mean value is entirely consistent 
with simply taking an arithmetic mean of metal- 
licities of either the Milky Way (-1.30 ±0.05 dex) 
or M31 (-1.14 ± 0.04 dex) globular cluster sys- 
tems. The dispersion about the mean is 0.88±0.10 
dex, larger than for the other two systems (~ 0.55 
dex). This is likely due to the sizeable uncertain- 
ties in the metallicities from Perelmuter, Brodie 
& Huchra (1995) and Brodie & Huchra (1991) 
rather than an intrinsically larger spread in the 
metallicities of the M81 globular clusters. 

3.2. Internal Reddening in M81 

Using the B— V color-metallicity relation estab- 
lished for Milky Way globular clusters by Couture 
et al. (1990), we computed metallicities predicted 
by their B — V colors for the clusters in the Keck 
sample. Object 12 was excluded from this anal- 
ysis (see Section 2). Figure 4 shows the photo- 
metric metallicities plotted against the metallic- 
ities derived from their spectra as described in 
Section 3.1. For reference we also show Milky 
Way clusters for which foreground reddening is low 
(E(B-v)< 0-4) and the line where clusters would 
lie if photometric and spectroscopic metallicities 
were in perfect agreement In the upper panel, the 
M81 cluster colors were corrected for foreground 
reddening only. Even after such a correction, the 
majority still lie above the line, i.e., their spec- 
troscopic metallicities are lower than their colors 
would predict. An obvious explanation is that the 
light from the M81 clusters suffers internal red- 
dening, which is likely to be variable across the 
disk of the galaxy. The range of color excess that 
would have to be applied to get all of the M81 clus- 
ters above the line to lie directly on it is 0.01 < 
E(B_v)< 0.24. Given both the intrinsic scatter and 
non-linearity of the true color-metallicity relation, 
we cannot know the precise reddening for the in- 
dividual clusters. However, the mean of the range 
quoted above (E(b_v)) =0.12 is likely a reasonable 
estimate of the average internal reddening in M81. 
This value is slightly lower than the value quoted 
by Perelmuter & Racine (1995) but to within the 



7 



errors they are consistent. The lower panel of Fig- 
ure 4 shows the results of applying an additional 
0.1 mag correction to the cluster colors and recom- 
puting the photometric metallicities. 

To investigate the degree to which reddening 
could affect our derived cluster metallicities, we 
further de-reddened the spectra of the three clus- 
ters which required the largest corrections to bring 
their photometric and spectroscopic metallicities 
into agreement and remeasured their metallici- 
ties. We found that de-reddening the spectra sys- 
tematically increases our derived mean metallic- 
ities, but the average change is only 0.05 dex, 
which is a fraction of the average statistical un- 
certainty of ~0.20 dex. The errors themselves 
also systematically increase, due to increased scat- 
ter in the metallicity values predicted by indi- 
vidual spectral features. This exercise demon- 
strates that spcictroscopic metallicity determina- 
tions are far less vulnerable to errors introduced 
by internal reddening, a significant problem for 
the method of determining metallicities by pho- 
tometry alone. Spectroscopy is especially advanta- 
geous when observing clusters in the disks of spiral 
galaxies where patchy dust can make the internal 
reddening highly variable. 

3.3. Individual Absorption Line Strengths 

Our spectra have a sufficiently high signal-to- 
noise ratio to compare the strengths of individual 
absorption features to those of Milky Way and 
M31 globular clusters. The complete suite of 
Lick/IDS indices (Gonzalez 1993; Trager 1998) 
for the globular clusters in our M81 sample are 
presented in Table 6a and Table 6b. Figure 5 
shows the strengths of G4300, NaD, TiOl, CN2, 
(Fe) and Mg2 plotted against [MgFe] strength 
for M81, Milky Way and M31 globular clusters. 
Gonzalez (1993) defines (Fe) as the arithmetic 
mean of the F( 333r) and Fe5270 indices, and 
[MgFe] as ^ (Fe) x Mgb. The relations between 
these indices and the [MgFe] index for the M81 
clusters in the Keck sample are generally the same 
as for Milky Way and M31 globular clusters. We 
had sufficient spectral coverage to measure the 
CN2 strengths for five M81 clusters, in which we 
see no evidence of a CN enhancement (see panel 
d of Figure 5). 

The Mg2 index appears to be slightly depressed 
in M81 clusters when compared to M31 and Milky 



Way clusters with similar [MgFe] values (sec panel 
e of Figure 5). In Figure 6 the values of Fe5270, 
Fe5335 and Fe5406 are plotted against Mgl, Mg2 
and Mgb for M81, M31 and Milky Way clusters. 
Mg2 appears depressed in M81's globular clusters, 
relative to Fe5335 but it is completely in line for 
the other two iron indices. The depression of Mg2 
in Figure 5 can be better described as an excess 
of Fe5335, causing enhanced [MgFe] values. But 
is it a real effect? We do not think it is likely. 
The strengths of the Fe5335 and Fe5270 indices, 
whose average makes up (Fe) and thus affect the 
value of [MgFe], can be influenced by contami- 
nation from other elements. Of all the Lick/IDS 
iron indices, only Fe5406 is a pure measure of iron 
(Trager 1998), and the Fe5406 index values of the 
M81 clusters are quite consistent with the Milky 
Way and M31 clusters for all three of the magne- 
sium indices. The relative abundances of magne- 
sium and iron in M81's globular clusters appear 
to be completely consistent with Milky Way and 
M31 globular clusters, at least at the upper end of 
the metallicity range. 

3.4. Ages 

In Figure 7 we plot the age-sensitive index 
H/3 against Mg2 for M81 globular clusters. For 
reference we also show these index values for 
individual Milky Way and M31 globular clus- 
ters, taken from the work of Trager (1998). 
Isochrones from the evolutionary-synthesis models 
of Worthey (1994) and from Fritze-v.Alvensleben 
& Burkert (1995) are overplotted. Although there 
are significant differences between the models (the 
Worthey models give lower H/3 absorption val- 
ues at a given age and metallicity, probably be- 
cause these models do not include blue horizontal- 
branch stars in their population-synthesis scheme) 
it is clear that all but one of the M81 clusters 
in our sample are old, i.e., their ages are simi- 
lar to those of Milky Way and M31 globular clus- 
ters. The exception is object 15, whose H/3 index 
value is quite high when compared to the other 
clusters in the sample. Strong H/3 absorption in 
globular clusters is often interpreted as indicat- 
ing low cluster metallicity, since the horizontal 
branches of metal-poor clusters are typically bluer 
than those of metal-rich clusters. However, the 
model tracks are isochrones, showing the behav- 
ior of the H/3 equivalent width for simple stellar 



8 



populations of a single ago, for a range of mctal- 
licities. The strong H/3 absorption in this cluster 
cannot be explained by metallicity alone. In ad- 
dition, the color of this cluster is much bluer than 
would be expected on the basis of its mean metal- 
licity. This object may in fact be a young (3 Gyr) 
globular cluster. Similar clusters have been ob- 
served in both M33 and M31 (Brodie & Huchra, 
1990; 1991). 

4. Cluster Kinematics 

Using the maximum likelihood method of Pryor 
& Meylan (1993) we find a mean line-of-sight ve- 
locity for the Keck sample of —105 ± 32kms~^ 
and a velocity dispersion of 123 ± 23kms~^. The 
mean velocity value differs significantly from both 
the heliocentric systemic velocity of M81 (—34 ± 
4kms~^) and the median cluster radial velocity 
found by Perelmuter, Brodie & Huchra (1995) of 
— 18 ± 32 kms"-^. We note, however, that the spa- 
tial extent of our sample is small and biased to 
one side of the galaxy (see Figure 1). The Keck 
sample shows no significant evidence for rotation, 
but the spatial bias could be obscuring its pres- 
ence. To increase the sample's size and spatial 
extent, we include the velocities of 25 confirmed 
clusters from the study of Perelmuter, Brodie & 
Huchra (1995) as well as the unpublished ve- 
locities of eight additional objects from Brodie & 
Huchra (1991). As mentioned in Section 3.1, we 
have matched all eight of the Brodie & Huchra 
objects with members of the Perelmuter & Racine 
(1995) candidate Hst (see Table 5). Based on their 
colors and/or proper motions, two of them are 
likely foreground stars. The remaining six ob- 
jects are confirmed M81 globular clusters, three of 
which were also observed by Perelmuter, Brodie & 
Huchra (1995). For these three we have adopted 
the error-weighted mean of the velocities from the 
two studies. We combined the three data sets into 
a sample containing 44 M81 clusters, which we 
henceforth refer to as the "full sample" . 

The full sample has a mean heliocentric line- 
of-sight velocity of —34 ± 23kms^^, in excellent 
agreement with M81's systemic velocity. The line- 
of-sight velocity dispersion of the full sample is 
133 ± 19kms~^, generally consistent with the ve- 
locity dispersions of the globular cluster systems of 
the Milky Way {a ~ 100; Armandroff, 1989; Da- 



Costa & Armandroff, 1995) and M31 {a - 150; 
Huchra, 1993; Huchra, Brodie & Kent, 1991). The 
line-of-sight dispersion of M104's globular cluster 
system is significantly larger (cr ^ 260; Bridges et 
al. 1997), but that is not unusual given its higher 
luminosity. 

4.1. Rotation 

Figure 8 shows the positions of the M81 glob- 
ular clusters marked with circles at the location 
of each cluster. The size and thickness of the 
circles are indicative of the magnitude and direc- 
tion of the cluster radial velocities. Receding clus- 
ters {thick-lined circles) tend to be located in the 
northwest, approaching clusters {thin-lined cir- 
cles) tend to lie more to the southeast. This sug- 
gests that rotation is present in the M81 globular 
cluster system in a sense that is consistent with the 
disk (Adler & Westpfahl 1996). A more quantita- 
tive analysis is demonstrated in Figure 9. Panel a 
shows the velocities (relative to the galaxy) of the 
full sample of clusters as a function of position an- 
gle with respect to the galaxy center. We used the 
maximum-likelihood method of Pryor & Meylan 
(1993) to determine a rotation velocity and rota- 
tion axis for the sample. The analysis assumes 
that the system is rotating as a solid body. Data 
points taken from true solid-body rotators would 
follow the sinusoids in Figure 9 quite closely, al- 
though some scatter would necessarily be present, 
since clusters at the same position angle will have 
a variety of perpendicular distances from the axis 
of rotation. The rotation velocity quoted for each 
sub-sample in Figure 9, Figure 10 and in Table 7 
is the linear speed of a point in the solid body 
which the sinusoid represents taken at the mean 
projected galacto centric radius of the sub-sample. 
For the full sample we find a rotation velocity of 
85 ± 32kms-i about an axis that is 17±28 de- 
grees east of north. This axis, shown by the arrow 
in Figure 8, is roughly aligned with the galaxy's 
minor axis and also the rotation axis for M81's gas 
(Adler & Westpfahl 1996). For comparison, the 
same analysis of 182 M31 globular cluster veloci- 
ties yields a rotation amplitude of 55 ± 14kms~"'^ 
and a line-of-sight dispersion of 138 ± 8kms~^. 
Here too, the rotation axis derived for the entire 
sample of M31 clusters is roughly aligned with 
M31's minor axis. 

Given the division of the Milky Way and 
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M31 globular cluster systems into kinematic sub- 
systems which differ in their spatial and chemical 
characteristics, it is of interest to look for evidence 
of such sub-populations in the globular cluster sys- 
tem of M81. Since the spectroscopic metallicities 
of about two- thirds of the M81 full sample have 
large uncertainties, the clusters' B— R colors must 
sufSce for assigning them to one metallicity pop- 
ulation or another. In the Milky Way globular 
clusters, the break between the metal-rich and 
metal-poor clusters corresponds to B — R~ 1.3, so 
we adopt this value for separating the M81 clus- 
ters into red (mctal-rich) and blue (metal-poor) 
sub-samples containing 33 and 11 members, re- 
spectively. We also examined a group of the 18 
reddest clusters (B — R> 1.5). We realize that 
these sub-samples represent only rough metallic- 
ity groups and that, because of internal reddening 
(see Section 3.2), some clusters which belong in 
the low-metallicity group could be "scattered" 
into the high-metallicity group, i.e. the red sub- 
samples may contain clusters that are not truly 
metal-rich. 

The results of our kinematic analysis are sum- 
marized Table 7. Figure 9 illustrates that rota- 
tion is present in both red sub-samples at the 2 a 
level, while the blue sub-sample shows no evidence 
for rotation. Additional analyses of sub-samples 
based on spatial characteristics alone showed that 
it is clusters at intermediate projected galacto- 
centric radii that contribute most strongly to the 
global rotation signature (Figure 10). As Ta- 
ble 7 shows, the precise value of the rotation ve- 
locity for intermediate-range clusters varies de- 
pending on the choice of radius. We found the 
most significant detection in the 3-7 kpc sub- 
sample (157±37kms^^), but the other reasonable 
choices of radius limits (2-6 kpc and 4-8 kpc) 
yielded strong detections also. The presence of 
strong rotation in MBTs globular clusters at in- 
termediate projected galactocentric radii invites 
a comparison to the Milky Way thick-disk clus- 
ters noted by ArmandrofF (1989) which were 
further described by Burkert & Smith (1997) 
as a "ring" of rapidly rotating metal-rich Milky 
Way clusters in a 2 kpc wide annulus at interme- 
diate galactocentric radius. Most of the M81 clus- 
ters in the intermediate-radius sample are mem- 
bers of the red, mctal-rich sample, i.e., most have 
B — R > 1.3. Unfortunately our sample is not 



large enough to make a direct comparison between 
metal-rich and metal-poor clusters in specific ra- 
dius ranges, but at the very least our data are 
consistent with M81's mctal-rich clusters at inter- 
mediate projected radii being associated with a 
thick disk in M81. 

One may wonder why such a population was 
not detected in the globular cluster system of M31, 
given the much larger sample of clusters for which 
precise velocities and reliable spectroscopic metal- 
licities have been measured. The main result of 
Huchra, Brodie & Kent (1991)'s kinematic study 
was that the inner metal- rich clusters in M31 were 
the most rapidly rotating. The analysis that Ash- 
man & Bird (1993) performed on the Huchra, 
Brodie & Kent (1991) data, showed strong sta- 
tistical support for the division of M31's globu- 
lar cluster system into disk and halo components. 
The method employed by Huchra, Brodie & Kent 
(1991) was quite different than the one we use 
here. They used each cluster's projected distance 
along the galaxy major axis to bin the metal-rich 
and metal-poor clusters into spatial sub-samples. 
Then, for each metallicity group they simply com- 
puted the difference in the average heliocentric 
radial velocity for clusters on either side of the 
galaxy, lying beyond a given distance along the 
major axis. They estimated the rotation veloc- 
ity to be one-half of this value. They found the 
difference between the metal-rich and metal-poor 
cluster rotation velocities to be most pronounced 
at very small distances along the galaxy major axis 
(< 1.5 kpc) with less pronounced differences per- 
sisting out to distances of roughly 8 kpc. They did 
not examine sub-samples in specific ranges of pro- 
jected radii along the galaxy major axis, so their 
analysis did not address the question of whether 
the rotation of clusters at intermediate projected 
radii might show even stronger differences between 
the metal-rich and metal-poor groups. 

We performed the maximum-likelihood kine- 
matic analysis employed for the M81 clusters on 
166 M31 clusters for which spectroscopic metal- 
licities exist (Barmby et al. 2000). Table 8 
shows that, as in the Milky Way, the most signifi- 
cant difference between the rotation of the metal- 
rich and metal-poor clusters occurs at interme- 
diate projected galactocentric radii. For all rea- 
sonable choices of radii, the; difference between 
the metal-rich and metal-poor clusters is quite 
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striking. Admittedly, treating the rotation axis 
as a free parameter greatly increases the chance 
for detecting rotation in a small sample. As a 
case in point, our analysis yields a weak detec- 
tion of rotation with approximately the same am- 
plitude but with very different rotation axes for 
metal-rich and metal-poor clusters lying inside 2 
kpc. However, most of the metal-poor sub-samples 
have poorly-constrained rotation axes and have 
rotation- velocity-to- velocity-dispersion ratios that 
are consistent with their being supported primar- 
ily by thermal motion rather than by rotation. 
(The metal-rich M31 clusters inside 2 kpc also 
have a low indicating that they may be 

associated with M31's bulge.) The metal-rich 
sub-samples all have tlicnr rotation axes roughly 
aligned with the galaxy's minor axis and they are 
rotating in the same sense as M31's gas (Loinard, 
Allen & Lequeux 1995). Still, as a check on the 
extent to which treating the rotation axis as a free 
parameter could influence the value of the rotation 
velocity derived for the M31 sub-samples, we fixed 
the rotation axis to be aligned with M31's minor 
axis and repeated the analysis. To within the un- 
certainties, this procedure little or no change in 
the rotation velocities for all sub-samples except 
for the metal-poor clusters inside 2 kpc, which 
exhibit no rotation at all about the galaxy's mi- 
nor axis. Thus the main result of Huchra, Brodie 
& Kent (1991), who assumed the clusters were 
rotating about the galaxy minor axis, is recovered 
here. However, our analysis also reveals a poten- 
tial thick-disk population among M31's metal-rich 
globular clusters, analogous to those found in the 
Milky Way and for which we see strong evidence 
in the globular clusters of M81. 

While nature has not provided a sufficient num- 
ber of metal-rich clusters residing at intermedi- 
ate projected galactocentric radii to ever con- 
struct samples large enough to allay concerns 
about small-number statistics, the fact that all 
three spiral galaxies whose globular cluster sys- 
tems have been studied in reasonable spectro- 
scopic detail show evidence for similar populations 
offers encouragement that they are part of analo- 
gous stellar populations, which are possibly the 
thick disk components of these galaxies. 



4.2. Projected Mass 

The velocity distribution of the globular clus- 
ters can be used to give an estimate of the 
mass of the parent galaxy. To this end we fol- 
low the method of Heisler, Tremaine & Bah- 
call (1985) and use the projected mass estimator 

We adopt a value of a = 1.5, as recommended by 
Heisler, Tremaine & Bahcall (1985) and assume 
a value of /p = The latter corresponds to 
the case that the clusters are on isotropic orbits. 
If the clusters are on radial orbits, then the de- 
rived projected mass will be lower by a factor of 
2. Based on the full sample of 44 clusters, the 
projected mass of M81 enclosed within 21.3 kpc 
is 4.0 ± 0.8 X 10" M©. The error quoted is the 
dispersion in the distribution of projected masses 
computed from 1000 simulated data sets gener- 
ated using the bootstrap method. This is only an 
approximation of the statistical error in the value 
of Mproj ■ It does not account for any error intro- 
duced by incorrect assumptions about the nature 
of the cluster orbits. 

Tenjes, Hand & Einasto (1998) quote a total lu- 
minosity for M81 of (Lb)m81 = 1-8 ± 0.3 x W° Lg, 
Using the mass derived above we determine 
a lower limit for M81's mass-to- light ratio of 
(M/Lb)m81 > 22, in solar units. 

Figure 11 shows a rough mass profile for M81 
based on masses derived using clusters inside pro- 
jected radii of 4, 8, 15 and 22 kpc. The specific 
values are listed in Table 9. The enclosed mass 
increases at large radius, indicative of a dark mat- 
ter halo around M81. Evidence for a dark matter 
halo was also detected by Perelmuter, Brodie & 
Huchra (1995). 

5. Discussion &: Conclusions 

We have obtained spectra of 16 globular clus- 
ter candidates in M81 with sufficient signal-to- 
noise to determine velocities and line-strengths for 
a number of absorption features. We also have 
supplemented our velocity data with velocities of 
28 additional clusters from previous work, in order 
to investigate velocity structures. Our conclusions 
are as follows: 

1. All 16 candidates we observed are confirmed 
as bona fide M81 globular clusters. One ap- 
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pears to have been undergoing a transient 
event at the time of observation. 

2. Our spectroscopic sample is biased to the 
metal-rich end of the metalhcity range. Due 
to the spatial location of the clusters near 
the center of the galaxy, the Keek sample 
is likely to be dominated by bulge and/or 
thick-disk clusters. The former seems more 
likely, given the lack of a strong rotation sig- 
nature in the Keck sample alone, but we note 
that rotation could be obscured by spatial 
bias in such a small sample. Based on the 
spectroscopically-determined metallicities of 
clusters in our sample, combined with metal- 
hcity estimates from previous work, we find 
the metallicity of the entire M81 globular 
chister system to be similar to the metal- 
licities of the Milky Way and M31 systems. 

3. Nearly all the clusters' ages, based on 
H/3 absorption-line strengths, are consistent 
with ages of Milky Way and M31 clusters. 
However, one cluster in the sample may be 
as young as a few Gyrs. Similar young clus- 
ters have been observed in both M31 and 
M33. 

4. Correlations between absorption line indices 
which have been established for globular 
clusters in the Milky Way and M31 hold for 
the M81 clusters as well, at least at the up- 
per end of the metallicity-range where the 
Keck sample probes. 

5. Based on the projected mass estimator, the 
mass of M81 is similar to the masses of the 
Milky Way and M31. The projected mass 
profile of M81 is consistent with the presence 
of a dark matter halo. 

6. The M81 cluster system appears to be kine- 
matically similar to the Milky Way and M31 
systems in that the metal-rich clusters show 
evidence for significant rotation while the 
metal-poor clusters exhibit halo-like kine- 
matics. The clusters which are rotating are 
doing so in the same sense as M81's disk and 
the global-rotation signature in our sample is 
driven primarily by clusters with projected 
galactocentric radii between 4 and 8 kpc. 
This may indicate the presence of a thick 



disk in M81. A similar analysis of the M31 
globular cluster system suggests that such a 
population is present in that galaxy as well. 

It is interesting to note that the annulus of 

rapidly-rotating, metal-rich clusters resides 
at the same galactocentric distance in all 
three galaxies discussed here. One inter- 
pretation is that the spatial and kinematic 
properties of these globular clusters do not 
obey galaxy scaling laws, but then the extent 
to which various components of the galaxies 
themselves follow scaling laws is not yet on 
firm footing. The thin (optical) disk of M31 
is larger than that of the Milky Way (Walter- 
bos & Kennicutt 1988; Binney & Tremaine 
1987), and M31's metal-rich clusters extend 
farther than their counterparts in the Milky 
Way. Unfortunately the radial scale length 
of the thick disk component, with which the 
rotating metal-rich clusters are thought to 
be associated, is only partially constrained in 
the Milky Way (see Buser, Rong & Karaali 
(1999) and references therein), and it has 
not been measured directly for either M31 or 
M81. Also, recent estimates of the dynam- 
ical masses of M31 and the Milky Way us- 
ing satellite galaxies, distant globular clus- 
ters and planetary nebulae as tracers sug- 
gest that, in spite of its larger appearance, 
M31 may actually be less massive than the 
Milky Way (Evans & Wilkinson 2000). Ob- 
viously more kinematic and metallicity data 
are needed for globular clusters in a vari- 
ety of spirals to determine if an annulus 
of rapidly rotating metal-rich globular clus- 
ters are present in all spirals and, if not, 
whether its presence depends on any intrin- 
sic or environmental property of the parent 
galaxy. 

In summary, the globular cluster system of M81 
is both kinematically and chemically very similar 
to the systems of the Milky Way and M31. This 
result lends support to the idea that universal pro- 
cesses govern the formation of globular cluster sys- 
tems around spirals. It may also be relevant to the 
formation of early-type galaxies. As has been dis- 
cussed extensively (Forbes, Brodie & Grillmaire 
1997; Ashman & Zepf 1998), it is difficuh to ex- 
plain the globular cluster systems of high specific 
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frequency ellipticals by a simple spiral-merger sce- 
nario, but this mechanism could be relevant in pro- 
ducing low specific frequency ellipticals. The blue 
sub-populations of globular clusters in such galax- 
ies, having originated in the spiral galaxy progen- 
itors, should be very similar. 
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Table 1 

Observed Cluster Candidates 











Aa 


Ad 


Up 








Slitlet 


ID 


(tC 720001 


(5f T2000') 




( PiTCVr\\T\l 




V 


B-V 


V-R 








(A A 






K') 


(s^ 
K^) 




(lOj 


1 


50285 


09:55:55.090 


69:00:56.14 


1.908 


-2.970 


3.530 


18.54 


1.10 


0.66 


2 


50304 


09:55:49.326 


69:01:15.29 


1.391 


-2.651 


2.994 


18.97 


1.02 


0.60 


3 


50359 


09:55:37.461 


69:02:07.58 


0.329 


-1.779 


1.810 


18.35 


1.04 


0.59 


4 


50378 


09:55:57.872 


69:02:23.07 


2.154 


-1.521 


2.637 


19.24 


1.05 


0.62 


5 


50418 


nq-55-54 642 


6q-02-52 57 


1 865 


— 1 030 


2 130 


18 45 


1 04 


57 


6 


50460 


09:55:51.466 


69:03:23.60 


L580 


-0^512 


L661 


18^80 


0^97 


0^46 


7 


50514 


09:55:48.157 


69:03:52.03 


1.284 


-0.039 


1.285 


19.05 


1.00 


0.54 


8 


50667 


09:55:22.295 


69:05:19.16 


-1.025 


1.414 


1.746 


18.01 


1.04 


0.56 


9 


50690 


09:55:21.587 


69:05:31.98 


-1.088 


1.627 


1.958 


18.76 


0.98 


0.48 


10 


50738 


09:55:30.277 


69:06:06.17 


-0.312 


2.197 


2.219 


20.08 


0.84 


0.68 


11 


50759 


09:55:35.841 


69:06:25.53 


0.184 


2.520 


2.526 


18.59 


0.82 


0.40 


12 


50782 


09 55 33.042 


69 06 39.88 


-0.066 


2.759 


2.760 


18.70 


1.00 


0.52 


13 


50787 


09:56:05.569 


69:06:42.97 


2.834 


2.810 


3.991 


19.12 


0.87 


0.57 


14 


50834 


09:55:25.403 


69:07:14.72 


-0.747 


3.340 


3.422 


19.03 


0.92 


0.51 


15 


50867 


09:55:51.995 


69:07:39.32 


1.622 


3.749 


4.085 


19.83 


0.75 


0.30 


16 


50889 


09:55:40.194 


69:07:30.82 


0.571 


3.608 


3.653 


18.73 


0.99 


0.62 


M81 Center 




09:55:33.780 


69:03:54.34 


0.000 


0.000 


0.000 









Note. — IDs, coordinates, magnitudes and colors were obtained from Perelmuter & Racine (1995). 

Slitlet number on LRIS mask. 
(^^ID number. 

Equatorial coordinates, J2000. 
(^'^^Position offsets from M81 center (arcminutes). 
(^)Projected radius from M81 center (arcminutes). 
Apparent V magnitude. 

Colors, uncorrected for reddening. 



Table 2 

Velocities for 16 M81 Globular Clusters 



Slitlet 
(1) 


ID 

(2) 


(kms 
(3) 




1 


50285 


-199 ( 


'33) 


2 


50304 


19 1 


27) 


3 


50359 


-131 1 


17) 


4 


50378 


-318 1 


31) 


5 


50418 


-63 ( 


33) 


6 


50460 


59 1 


32 


7 


50514 


-335 ( 


28) 


8 


50667 


-144 ( 




9 


50690 


-261 ( 




10 


50738 


22 1 


'39 


11 


50759 


11 1 


24) 


12 


50782 


-200 1 


17) 


13 


50787 


-87 ( 


^49) 


14 


50834 


69 1 




15 


50867 


-13 1 




16 


50889 


-100 (25) 



Slitlet number on LRIS 
mask. 

(^^ID number from Perel- 
muter & Racine (1995). 

^^■'Heliocentric line-of-sight 
velocity (error in parenthesis). 



Table 3 

BRODIE & HUCHRA METALLICITY 



— nm — 
but 


T 

A 


TTlTfTj 

CNB 


CNR 


7n 

G 


TTT — TT 

MgH 


TTT — o 

Mg2 


Fe52 


Nal 


— TTn — 7TTTT 

([le/HJ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 

J. 












1 iq4 






-1 OOT 












0.073 


0.038 


0.105 


0.064 


0.369 


9 












n 81 9 

-u.o±z 




-U. ( ( v 


n 7n7 




















n ifi7 


Q 
O 








u.uou 


-u.oio 


-u.ziu 


n 91 R 


n 4(17 


n 91 1 










0.191 


067 


036 


100 


063 


193 


4 










-0.470 


-0.311 


-0.505 


-0.887 


-0.407 












0.122 


0.065 


0.177 


0.111 


0.088 


5 










-1.210 


-1.062 


-0.972 


-0.821 


-1.086 












0.072 


0.038 


0.108 


0.067 


0.091 


6 










-1.785 


-1.415 


-1.262 


-1.383 


-1.493 












0.099 


0.052 


0.151 


0.095 


0.206 


7 










-1.057 


-0.971 


-0.794 


-1.151 


-0.955 












0.125 


0.066 


0.187 


0.116 


0.098 


8 


0.104 


-0.169 


-0.706 


-0.615 


-0.710 


-0.744 


-0.754 


-1.082 


-0.698 




0.057 


0.350 


0.070 


0.122 


0.053 


0.028 


0.079 


0.047 


0.058 


9 


-0.233 


-0.717 


-1.125 


-1.108 


-1.254 


-1.375 


-1.029 


-1.488 


-1.212 




0.133 


1.062 


0.127 


0.213 


0.092 


0.048 


0.136 


0.082 


0.133 


10 






-1.399 


-1.811 


-1.618 


-1.209 


-0.750 


-1.286 


-1.322 








0.758 


0.928 


0.348 


0.183 


0.527 


0.380 


0.356 


11 






-2.222 


-0.558 


-1.691 


-1.193 


-1.193 


-1.745 


-1.114 








0.092 


0.175 


0.069 


0.036 


0.104 


0.070 


0.409 


13 










-1.074 


-1.099 


-0.948 


-1.568 


-1.055 












0.149 


0.078 


0.208 


0.135 


0.062 


14 


-0.210 


-0.315 


-1.140 


-1.028 


-1.179 


-1.171 


-1.023 


-1.385 


-1.107 




0.244 


2.468 


0.142 


0.239 


0.096 


0.051 


0.145 


0.094 


0.074 


15 








-0.080 


-1.981 


-1.251 


-0.400 


-1.066 


-1.014 










0.468 


0.173 


0.092 


0.270 


0.191 


0.713 


16 










-0.687 


-0.630 


-0.739 


-0.804 


-0.674 












0.100 


0.053 


0.145 


0.091 


0.044 



Note. — There are two lines for each cluster. The first line contains the metallicity 
measurements. The second line contains the uncertainties in the metallicity measure- 
ments. The uncertainties for the individual metallicity measurements were determined 
by determining a "total error" for the index measurements, (equation 19 of Brodie & 
Huchra, 1990) and then propagating that total error through their derived linear index- 
metallicity relationships. The uncertainty quoted for the weighted mean metallicity for 
each cluster is weighted dispersion about weighted mean, where the weights were derived 
using equation 20 of Brodie & Huchra, 1990. 

^^^Slitlet number on LRIS mask. 

(^"^^Metallicities predicted by various absorption line strengths, as described in Brodie 
& Huchra (1990). 

(lO)Error -weighted mean, does not include A or Nal. 



Table 4 

Indices Used to Predict Metallicity 



yiitlet 


A 


CiNii 


(JiNR 


(J 


MgH 


Mg2 


l''e52 


J\al 


(1) 




X'-'J 


(A) 


(5) 


(6) 




(8) 


(9) 


1 










009 


103 


070 


139 












0.004 


0.004 


0.005 


0.004 


2 










0.067 


0.141 


0.062 


0.110 












0.005 


0.005 


0.007 


0.007 


3 








0.218 


0.064 


0.201 


0.092 


0.135 










0.017 


0.003 


0.004 


0.005 


0.004 


4 










0.067 


0.192 


0.078 


0.102 












0.006 


0.007 


0.009 


0.008 


5 










0.031 


0.116 


0.055 


0.106 












0.004 


0.004 


0.005 


0.005 


6 










0.003 


0.080 


0.040 


0.068 












0.005 


0.005 


0.007 


0.007 


7 










0.038 


0.125 


0.063 


0.084 












0.006 


0.007 


0.009 


0.008 


8 


0.796 


0.269 


0.066 


0.161 


0.055 


0.148 


0.065 


0.088 




0.018 


0.053 


0.010 


0.011 


0.003 


0.003 


0.004 


0.003 


9 


0.690 


0.186 


0.009 


0.118 


0.028 


0.084 


0.052 


0.060 




0.042 


0.162 


0.017 


0.019 


0.004 


0.005 


0.007 


0.006 


10 






-0.028 


0.056 


0.011 


0.101 


0.066 


0.074 








0.103 


0.081 


0.017 


0.018 


0.026 


0.026 


11 






-0.140 


0.166 


0.007 


0.103 


0.044 


0.043 








0.013 


0.015 


0.003 


0.004 


0.005 


0.005 


13 










0.037 


0.112 


0.056 


0.055 












0.007 


0.008 


0.010 


0.009 


14 


0.698 


0.247 


0.007 


0.125 


0.032 


0.105 


0.052 


0.067 




0.077 


0.375 


0.019 


0.021 


0.005 


0.005 


0.007 


0.006 


tl5 








0.208 


-0.007 


0.097 


0.083 


0.090 










0.041 


0.008 


0.009 


0.013 


0.013 


tio 










0.056 


0.159 


0.066 


0.108 












I). 005 


0.1)05 


0.007 


0.006 



Note. — There are two lines for each cluster. The first line contains the 
index measurements. The second line contains the uncertainty (photon er- 
ror) in the index measurements. Col.(l): Slitlet number on LRIS mask. 
Cols.(2)-(9): Absorption-line index values, measured as prescribed in Brodie 
& Huchra 1990. Uncertainties measured as described in text. 



Table 5 

Data for Additional M81 Clusters 













LAO 


Hp 








Vlos 










ID 


a(J2000) 


(5(J2000) 


(arcmin) 


(arcmin) 


(arcmin) 


V 


B- V 


V-R 


(kms' 




[Fe/H] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


R9 = 


HS18 = 50785** 


9:54:35.88 


69:06:43.1 


-5.161 


2.812 


5.878 


19.08 


0.86 


0.56 


229 1 


'66; 


1 -1.00 ( 


'0.84; 




HSl = 


= 40312 


9:54:52.01 


69:19:45.3 


-3.686 


15.849 


16.272 


19.11 


0.81 


0.55 


258 1 


27 


1 -2.10 ( 


0.97 




HS6 = 


= 60072 


9:55:08.08 


68:49:56.5 


-2.320 


-13.964 


14.156 


19.61 


0.80 


0.48 


122 1 


[77] 


1 -1.77 ( 


0.83, 




HS2 = 


= 50864* 


9:56:25.74 


69:07:27.8 


4.629 


3.558 


5.838 


19.62 


0.26 


0.26 


180 1 




\ -2.09 ( 




1 


R12 = 


= 60126* 


9:55:26.97 


68:44:53.6 


-0.617 


-19.013 


19.023 


19.16 


0.99 


0.59 


-37 1 




1 0.19 1 






R14 = 


= HS26 = 60012 


9:55:42.09 


68:55:00.6 


0.747 


-8.896 


8.927 


20.01 


0.97 


0.62 


-129 1 


;57; 


1 0.40 ( 


'2.65; 




HS35 


= 50960 


9:55:52.06 


69:08:18.7 


1.627 


4.406 


4.697 


18.49 


0.86 


0.54 


180 1 


59 


1 -1.43 1 


0.40 




HS19 


= 50415** 


9:56:20.58 


69:02:49.0 


4.184 


-1.089 


4.323 


19.24 


0.85 


0.48 


-218 1 


^51, 


1 -1.46 1 


^0.52; 





*Likcly to be a Galactic star, rejected from the sample. 

**Also observed by Perelmuter et. al, 1995 (See Section 3.1. 

(i)lD from Perelmuter & Racine (1995). 

(2'^)Equatorial Coordinates (J2000). 

(^'^^ Position offsets from M81 center (arcminutes) 

(^)Projected radius from M81 center (arcminutes). 

Apparent V magnitude. 
(8.9) Colors, uncorrected for reddening. 

Heliocentric, line-of-sight velocity (J. Huchra, private communication). Errors in parentheses. 
(^^'Metallicity from Brodie & Huchra (1991). 



Table 6a 
Lick IDS Index Values 





CNl 
{^) 


CN2 
(.3J 








ra4455 

\l) 


i''c4531 




(10) 




l^ 9^ 


M ffi 
(13) 


1 


















1.967 


3.500 


0.011 


0.10 




















0.197 


0.414 


0.004 


0.004 


2 
















0.790 


1.903 


3.688 


0.068 


0.139 


















0.950 


0.268 


0.578 


0.005 


0.005 


3 








6.651 


3.510 


2.012 


6.974 


4.097 


1.929 


5.226 


0.066 


0.199 










0.419 


0.600 


0.267 


0.374 


0.576 


0.175 


0.372 


0.003 


0.004 


4 


















2.246 


5.447 


0.069 


0.188 




















0.318 


0.669 


0.006 


0.006 


5 


















1.868 


3.357 


0.031 


0.114 




















0.187 


0.407 


0.003 


0.004 


6 
















3.083 


2.278 


2.166 


0.002 


0.080 


















0.852 


0.249 


0.563 


0.005 


0.005 


7 














4.119 


3.279 


1.832 


2.684 


0.040 


0.123 
















0.700 


1.058 


0.324 


0.726 


0.006 


0.007 


8 


0.068 


0.105 


1.090 


4.955 


3.919 


1.281 


2.942 


2.128 


1.898 


3.986 


0.055 


0.146 




0.009 


0.012 


0.171 


0.284 


0.411 


0.188 


0.294 


0.446 


0.135 


0.292 


0.003 


0.003 


9 


0.013 


0.057 


0.256 


3.675 


1.721 


1.102 


2.465 


1.040 


1.807 


3.215 


0.029 


0.084 




fi m 7 




U.OX X 






U.OO I 




n SI n 


n 9^7 


U.fJ 




U.UU(J 


10 


-0.038 


-0.015 


0.668 


1.484 


4.887 


1.275 


1.956 


2.117 


1.769 


6.273 


0.010 


0.100 




0.102 


0.127 


1.623 


2.420 


2.811 


1.262 


1.973 


2.949 


0.896 


1.842 


0.017 


0.018 


11 


-0.135 


-0.103 


1.180 


5.530 


4.098 


1.604 


2.696 


2.308 


1.766 


2.691 


0.007 


0.097 




0.012 


0.016 


0.211 


0.404 


0.539 


0.237 


0.391 


0.579 


0.177 


0.394 


0.003 


0.004 


13 




















2.946 


0.036 


0.112 






















0.850 


0.007 


0.008 


14 


0.004 


0.035 


0.544 


3.773 


1.394 


0.834 


2.657 


1.540 


2.167 


3.116 


0.032 


0.104 




0.019 


0.026 


0.343 


0.582 


0.850 


0.363 


0.552 


0.845 


0.250 


0.547 


0.005 


0.005 


15 






1.178 


6.072 


7.390 


1.291 


4.693 


-0.874 


4.509 


6.897 


-0.005 


0.099 








0.649 


1.050 


1.457 


0.711 


1.051 


1.544 


0.407 


0.888 


0.008 


0.009 


16 














3.681 


5.688 


1.643 


4.092 


0.057 


0.157 
















0.689 


0.891 


0.270 


0.560 


0.005 


0.005 



*^^^Slitlet number on LRIS mask. 

^^"^^^Absorbption-line index values measured as described in Trager 1998 and references therein. 



Table 6b 
Lick IDS Index Values, continued 



Elit Mgb hhb270 Fe5335 Fe5406 Fe5709 Fc5782 NiD TTUl H3F H^F HM mX 



(1) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


X 


9 Q3fi 


2.529 


2 063 

Zj. WDO 


1 269 


785 


762 


3 765 

O. I \jO 


021 


056 














1 73 

VJ. J. 1 (J 


217 

VJ.Zi J. 1 


149 


097 


094 


113 


003 

\J.\J\JO 


002 










2 




2 1 59 


2.io4 


1 026 


783 


547 


3 026 


032 


027 














246 


305 


21 3 

U<Zi J.O 


140 


1 37 


1 71 

U. -L 1 -L 


004 

\J.\J\J'-t 


003 












S Q27 


3 351 


3 059 


1 633 


941 


0.822 


3 604 


037 


078 




-1.112 




-4 977 






161 


201 


0.l42 


095 


092 


0.112 


003 

\J.\J\JO 


002 




255 




462 


4 


4.224 


2 901 


2 077 


1 675 


841 


969 


2 893 


025 


060 












0.27A 


0.289 


0.356 


0.245 


0.167 


0.160 


0.203 


0.005 


0.004 












2.444 


1 903 




754 


640 


682 


2 933 

Zj .iJOO 


019 


028 














180 


0.226 


1 59 


104 


100 


0.122 


003 

\J .yJyJO 


002 














1 398 


1 507 

-L .'J\J t 


\J .<Ji'J'.J 


3^1 


267 

\J . Z(U ( 


1 890 


012 


004 














2^5 


323 


0.224 


149 




180 


004 


003 










7 

t 


2 737 


2.233 


2 435 


L441 


770 


513 

\J .O J^O 


2.326 


037 


046 












n 293 


31 2 


387 


0268 


1 78 


1 74 


f) 21 ^1 


005 


004 










8 


2.635 


2.387 


2.081 


1.133 


0.780 
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*^^^Slitlet number on LRIS mask. 

^^^~^^^Absorbption-line index values measured as described in Trager 1998 and references therein. 



Table 7 

Kinematics of M81 Globular Clusters 











Color 
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Vrot 




N 


T? T-> Ran crc» 








(kms" 


(kms" 




(kms 






(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


u 


44 


Full Range 


7.9 




All 


-34 (23^ 


133 (19^ 


85 (32: 


17 (28^ 




33 


Full Range 


5.8 


B 


-R> 1.3 


-48 




143 


?^} 


93 




18 (33; 




11 


Full Range 


12.4 


B 


-R< 1.3 


-11 




11 




21 




143 (52^ 




18 


Full Range 


5.2 


B 


-R> 1.5 


-104 




121 




128 




-11 39; 




16 


Rpgc < 4kpc 


2.8 




All 


-103 




135 




89 




-17 (35; 




22 


2 < Rpgc < 6 kpc 


3.6 




All 


-56 (33; 


120 (25; 


119 (52; 


14 (28; 




17 


3 < Rpgc < 7 kpc 


4.2 




All 


-11 (38; 


103 (3o; 


157 (37; 


27 (20; 




13 


4 < Rpgc < 8 kpc 


5.3 




All 


1 (46: 


109 (37^ 


168 (42: 


45 (25; 




15 


Rpgc > 8 kpc 


14.2 




All 


13 (34; 


125 (3o; 


49 (52; 


31 (88; 





(^^Number of clusters in sample. 
Spatial selection criteria for sample. 

Average Rp of sample (assumes distance to M81 of 3.7 Mpc). 

Color selection criteria for sample. 
(^^Mean heliocentric line-of-sight velocity of sample (errors in parentheses). 

Velocity dispersion, corrected for rotation (errors in parentheses). 

Amplitude of sample rotation (errors in parenthesis). 
'^^•'Best fit position angle of rotation axis, positive angle is east of north. Errors in parentheses. 



Table 8 

Kinematics of M31 Globular Clusters 



Spatial 


Vrot 




t'rot 




Group 


(kms-i) 


(kms ^) 


n 


N 


(1) 


(2) 


(3) 


(4) 


(5) 



Rpgc < 2 kpc 
MR 
MP 

2 < Rpgc < 6 kpc 
MR 
MP 

3 < Rp < 7kpc 

MR 
MP 

4 < Rp < 8 kpc 

MR 
MP 

Rp < 8 kpc 
MR 
MP 



96 (56) 
81 (60) 

127 (40 
20 (28 

143 (48 
37 (33 

144 (68 
65 (34 

49 (39 
72 (26 



152 (36) 

153 (68) 

122 (40 
134 (15 

119 (43 
144 (16 

129 (50 
139 (21 

56 (27 
114 (13 



109 (56) 
-113 (51) 

119 (24) 
164 (108) 

128 (22 
136 (54 

87 (39 
121 (31 

96 (44) 
98 (29) 



17 
15 

16 
45 

14 
43 

11 
39 

12 
36 



(^^Limits of annulus defining sample. Assumes distance of 744 
kpc to M31. 

(^^Best fit rotation velocity, errors in parentheses 

^^■^ Velocity dispersion (corrected for rotation), errors in paren- 
theses. 

(■^^Best fit rotation axis (positive angle is east of north), errors 
in parenthesis 

Number of clusters in sample. 

Note. — MR=Metal-rich clusters, those with [Fe/H] >= 
—0.8. MP indicates clusters with metallicity less than this value. 



Table 9 
Projected Mass of M81 



N (kpc) (IQi^ Mq) 
(1) (2) (3) 



15 


3.9 


1.2 ( 


;o.8) 


30 


7.8 


2.4 ( 


;o.7) 


38 


14.7 


3.1 ( 


;o.7) 


46 


21.3 


4.0 1 


[0.8) 



(""^^ Number of points in 
sample. 

Maximum projected 
galactocentric radius in 
sample. 

^^■^ Projected mass interior 

to Rmax* 

Statistical uncertainty 
in projected mass. Error 
(in parentheses) does not 
include errors potentially 
introduced by assumptions 
about the cluster orbits. 
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Fig. 1. — Positions of all spectroscopically-confirmcd M81 globular clusters, relative to the center of M81. 
The ellipse shows the location and orientation of M81's disk, taken from Perelmuter, Brodie & Huchra (1995) 
Solid circles are members of the Keck sample. Open circles are clusters studied by Perelmuter, Brodie & 
Huchra (1995), and open squares are clusters from Brodie & Huchra (1991) for which velocities have been 
measured (J. Huchra, private communication.) Open triangles indicate clusters studied by both Brodie & 
Huchra (1991) and Perelmuter, Brodie & Huchra (1995). 
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Fig. 2. — Spectra of 16 M81 globular cluster candidates, obtained with 4500 seconds of integration with 
LRIS on the Keck I telescope. The spectra have been corrected for foreground reddening and shifted to 
the rest frame. Continuum shapes are reliable only between 3900A and 6500A. All are confirmed as M81 
globular clusters, although object 12 appears to have been observed during a transient event in the cluster 
(see Section 2). 
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Fig. 3. — The B — R color distribution of the Keck sample (upper panel) and the sample of Perelmuter, 
Brodie & Huchra (1995) (lower panel) compared to the color distribution of the entire list of M81 cluster 
candidates from Perelmuter & Racine (1995). The Keck sample is clearly biased to the redder, metal-rich 
end of the spectrum. The Perelmuter, Brodie & Huchra (1995) sample has a more uniform color distribution 
(see Section 3.2). 
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Fig. 4. — Photometric vs. spectroscopic metallicities for M81 (filled circles) and Milky Way (open circles) 
globular clusters. Photometric metallicities were computed from B — V colors using the color- metallicity 
derived in Couture et al. (1990). In the upper panel, the colors were corrected for foreground reddening 
only, which was estimated to be ~ 0.1 (Perelmuter & Racine 1995). The lower panel shows that de-reddening 
the M81 cluster colors by an additional 0.1 mag brings the photometric metallicities of the M81 clusters in line 
with those of the Milky Way (but see Section 3.2). Only Milky Way clusters for which foreground reddening 
is low (E(B_v)< 0.4) are included. The reddening and metallicity values for the Milky Way clusters were 
obtained from Harris (1996). 
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Fig. 5. — Selected Lick IDS indices plotted against the metallicity-sensitive index [MgFe] for M81 (filled 
circles), M31 (open squares) and Milky Way (open triangles) globular clusters. The error bars show the 
typical sizes of the statistical (photon) errors on the M81 cluster indices. The correlations between specific 
indices observed for Milky Way and M31 globular clusters are also present in the M81 globular cluster 
system. The apparent depression of Mg2 with respect to [MgFe] (panel e) is likely due to contamination of 
the Fe5335 index. See Section 3.1 and Figure 6 for more detail. 
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Fig. 6. — Iron indices plotted against magnesium indices for M81 (filled circles). M31 (open squares) and 
Milky Way (open triangles) globular clusters. The error bars show the typical sizes of the statistical errors 
for the M81 data. The Fe5335 index appears enhanced, but both Fe5335 and Fe5270 can be affected by the 
strength of absorption by other elements (Trager 1998). Fe5406 is a pure measure of the iron abundance, and 
based on this index we conclude that relative abundances of magnesium and iron in M81 globular clusters at 
the metal-rich end of the metallicity range is consistent with the abundances for globular clusters in M31 
and the Milky Way. The apparent depression of Mg2 in panel e of Figure 5 is likely due to contamination 
artificially enhancing the Fe5335 index. 




Fig. 7. — The age-sensitive index H/3 plotted against the metallicity-sensitive index Mg2 for M81 (filled 
circles), M31 (open squares) and Milky Way (open triangles) globular clusters. Overplotted as solid lines 
are Worthey (1994) isochrones of ages 1.5 and 3 Gyrs (spanning [Fe/H] -0.225 to 0.5) and 8 and 17 Gyrs 
(spanning [Fe/H] of -2.0 to 0.5). The dashed lines are isochrones of ages 1.5, 3, 8 and 15 Gyrs from the 
models of Fritze-v.Alvensleben & Burkert (1995), spanning metal-content of [Fe/H] = —1.3 to solar. With 
the exception of object 15, the M81 cluster ages are consistent with the ages of M31 and Milky Way globular 
clusters, i.e., they are old. 
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Fig. 8. — Velocities and positions of the full sample of M81 globular clusters. Thick-lined circles are receding, 
thin-lined circles are approaching. The ellipse shows the location of the MSl's disk, taken from Perelmuter, 
Brodie & Huchra (1995). The qualitative impression is that the M81 cluster system is rotating. The arrow 
shows the rotation axis for the full sample (see Section 4.1). 
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Fig. 9. — Quantitative assessment of rotation in M81 globular cluster sub-samples based on color. Panel a 
shows rotation for the full range of colors. The rotation velocities were derived using the method of Pryor 
& Meylan (1993). The vertical axis is the cluster radial velocity relative to the mean velocity of the full 
sample. The horizontal axis is the position angle between a cluster and the galaxy center, where a PA of 
degrees is directly north-south and positive PA is cast of north. The best fitting sinusoids to these data 
represent solid-body rotation and place a lower limit on the rotation velocity of each sample. Rotation is 
strongly suggested in both of the red, metal-rich sub-samples. The reddest sub-sample also has the lowest 
velocity dispersion (see Table 7). The blue, metal-poor sub-sample shows no significant evidence for rotation. 
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Fig. 10. — Quantitative assessment of rotation in M81 globular cluster sub-samples based on projected 
galactoccntric radius. Only the sub-sample of clusters with intermediate projected galactocentric radii show 
evidence for rotation. See Section 4.1. 
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Fig. 11. Projected mass vs. projected galactocentric radius for M81, derived using radial velocities of M81 
globular clusters (see Section 4.2). The increase in mass at large radius is consistent with the presence of a 
dark-matter halo. 



